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An elementary quantum network operation involves storing a qubit state in an atomic quantum 

in"; 

memory node, and then retrieving and transporting the information through a single photon excita- 

o . 

tion to a remote quantum memory node for further storage or analysis. Implementations of quantum 

>: 

network operations are thus conditioned on the ability to realize such matter-to-light and/or light-to- 
matter quantum state mappings. Here, we report generation, transmission, storage and retrieval of 



■ single quanta using two remote atomic ensembles. A single photon is generated from a cold atomic 

>: 

I" . ensemble at Site A via the protocol of Duan, Lukin, Cirac, and Zoller (DLCZ) [1] and is directed to Site 
t-H ■ 

B through a 100 meter long optical fiber. The photon is converted into a single collective excitation via 



ly-^ the dark-state polariton approach of Fleischhauer and Lukin [2j. After a programmable storage time 

o ; 

the atomic excitation is converted back into a single photon. This is demonstrated experimentally, for 

Oh' 

a storage time of 500 nanoseconds, by measurement of an anticorrelation parameter a. Storage times 

^— > ' 

exceeding ten microseconds are observed by intensity cross-correlation measurements. The length of 
■ the storage period is two orders of magnitude longer than the time to achieve conversion between pho- 

> : 

•rH . tonic and atomic quanta. The controlled transfer of single quanta between remote quantum memories 

X; 

^ . constitutes an important step towards distributed quantum networks. 

_ _ i 

A quantum network, consisting of quantum nodes and interconnecting channels, is an outstanding goal of quantum 
information science. Such a network could be used for distributed computing or for the secure sharing of information 



between spatially remote parties 



While it is natural that the network's fixed nodes (quantum memory 
elements) could be implemented by using matter in the form of individual atoms or atomic ensembles, it is equally 
natural that light fields be used as carriers of quantum information (flying qubits) using optical fiber interconnects. 
The matter-light interface seems inevitable since the local storage capability of ground state atomic matter cannot 
be easily recreated with light fields. Interfacing material quanta and single photons is therefore a basic primitive of a 
quantum network. 



The potential of atomic ensembles to serve as quantum memories has recently attracted considerable attention 

notfnnn 

[U, 12, La 121 UH 1111 , spawning two distinct lines of research. In one, using the physics of "slow light" propagation 
in an opticall y th ick atomic ensemble, weak coherent laser pulses have been stopped and retrieved in a controlled 
fashion 3 3 3- In the other, motivated by the seminal proposal of Duan, Lukin, Cirac, and Zoller (DLCZ) 
correlated pairs of photons and single photons have been produced from an atomic ensemble 3 3 17 , 3 3 20 1 . 
Collective atomic qubits, atom-photon entanglement, and quantum state transfer from atomic to photonic qubits have 
also been demonstrated . These initial experimental demonstrations within the DLCZ paradigm were beset by 
short atomic coherence times, of the order of the laser pulse length. In contrast, recent advances in atomic ensemble 
research 3 allow for long quantum memory times, in excess of ten microseconds in the present work, more than two 
orders of magnitude longer than the duration of the laser pulses involved in the protocols. 

Here we report the synthesis of these two lines of research by demonstrating the generation, transmission, storage 
and retrieval of single photons using remote atomic ensembles as quantum memories. The essential ingredient which 
enables the completion of this synthesis, and which we report here, is the ability to convert single photons into single 
collective atomic excitations. In our experiment the remote quantum memories are based on cold atomic clouds of 
85 Rb confined in magneto-optical traps (MOTs) at Sites A and B, as shown in Fig. 1. Sites A and B are physically 
located in adjacent laboratories, with a 100 meter long single-mode optical fiber serving as the quantum information 
channel. 

Our protocol begins with the generation of single photons at Site A, using an improved version of the DLCZ 
approach in the off- _„fi, s „ QQ X, fifi» c^e, tos t fie field * Sfie B 

where an optically thick atomic ensemble is prepared in level \b) (right inset in Fig. 1). The signal field propagation 
in the atomic medium is controlled by an additional laser field (control) through the process of electromagnetically- 
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induced transparency (EIT) j23ll'2-l| . As we deal with an unpolarized atomic ensemble, we must take into account the 
Zeeman degeneracy of the atomic levels. Choosing the same circular polarizations for both the probe and the control 
fields allows us to retain transparency, as discussed in more detail in the Supplementary Information. In Fig. 2 we 
show the EIT transmission spectrum recorded for a coherent laser probe field instead of the signal field. Evidently, 
in the absence of the control light the probe field is absorbed by the optically thick sample. With the addition of the 
cw control field, the medium is rendered transparent around the \b) «-> \c) transition resonance A = 0. 

The control field strongly modifies the group velocity of the signal field. For a time-dependent control field, a 
strong reduction of the group velocity of the propagating signal field can be understood in terms of a coupled matter- 
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FIG. 1: A schematic diagram of our experimental setup demonstrating generation, transmission, storage and 
retrieval of single photon excitations of the electromagnetic field. Two atomic ensembles at Sites A and B are connected 
by a single-mode fiber. The insets show the structure and the initial populations of atomic levels for the two ensembles. All the 
light fields responsible for trapping and cooling, as well as the quadrupole magnetic fields in both MOTs, are shut off during the 
period of the protocol. The ambient magnetic field at each Site is compensated by three pairs of Helmholtz coils (not shown) . 
Correlated signal and idler fields are generated at Site A. The signal field is transmitted via optical fiber from Site A to Site 
B, where it is converted to atomic excitation, stored for a duration T s , and subsequently retrieved. A Hanbury Brown- Twiss 
setup consisting of a beamsplitter BS and two detectors D2 and D3, together with detector Dl for the idler field, are used to 
verify the single photon character of the retrieved field. 



light field excitation known as a "dark-state polariton." By adiabatically switching off the control field, the coupled 
excitation can be converted into a pure atomic excitation, i.e., the signal field is "stopped" 2, wj]. An important 
condition to achieve storage is a sufficiently large optical thickness of the atomic sample, which enables strong spatial 
compression of the incident signal field In our experiment the measured optical thickness d ~ 8. Fig. 3 compares 
our observations with the predictions of a theoretical model. Fig. 3a compares the propagation of the signal pulse in 
vacuum and in the atomic medium under conditions of EIT with a cw control field. The observed pulse delay under 
conditions of EIT is about 20 ns, corresponding to more than three orders of magnitude reduction in group velocity. 
Fig. 3b shows the effect of turning off the control- storage field when the signal pulse is approximately centered in the 
medium, and the subsequent retrieval of the signal field when the control-retrieval field is switched back on after a 500 
ns storage time. Fig. 3c shows retrieval after a storage time of 15 /is. Qualitative agreement of the pulse shapes has 
been obtained in our theoretical analysis of the protocol using the full Zeeman structure of the atoms and a classical 
description of the signal field (Fig. 3d-f). 

In order to verify the single-photon character of the signal field (a) without storage, and (b) with storage and 
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FIG. 2: Measured transmission spectra of a coherent probe field as the function of probe detuning in the 
presence of, and absence of, EIT. Data are taken using 700 ns long coherent laser pulses. T is the intensity transmittance, 
A is the probe detuning and F is the decay rate of level |c). In the absence of control field (circles) the probe is strongly 
absorbed near resonance, whereas with the control field on (diamonds) the medium becomes transparent. Each probe pulse 
contains on average 0.3 photons. Each data point is an average of 2 x 10 5 experimental trials. The optical thickness d = 8 and 
the control field Rabi frequency Q — 3F are used to obtain the solid curves, based on the theoretical model discussed in the 
Supplementary Information. 

retrieval, we use a Hanbury Brown-Twiss detection scheme, employing a beamsplitter followed by two single photon 
counters, as shown in Fig. 1 |25| . To provide such characterization, we note that classical fields must satisfy a criterion 



a > 1 based on the Cauchy-Schwarz inequality 
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For an ideally prepared single photon state a — > 0. Here 
the anticorrelation parameter a is a function of the storage time T s , and is given by the ratio of various photoelectric 
detection probabilities which are measured by the set of detectors D1,D2 and D3 (described in the Methods section): 



,rp \ P1P123 n , 
a(T„) = . (1) 

P12P13 

As an auxiliary measure of signal-idler field correlations, and as a way to quantify the quantum memory storage 
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In 



time, we also evaluate the normalized intensity cross-correlation function g S i = (pi2 +Pi3)/[pi(p2 +P3)] 
particular, it serves to estimate the total efficiency and background levels in the experiment, since g S i is, by definition, 
independent of efficiencies whereas pi is proportional to the overall idler channel efficiency. 

First we measure g S i and a without storage at Site B (i.e., with no atomic sample in place), and the results are 
displayed in Fig. 4, a and b, respectively. Next we add an optically thick atomic sample at Site B, and perform 
storage of duration T s = 500 ns and subsequent retrieval of the signal field, with results shown in Fig. 4, c and d, 
respectively. No correction for background or dark counts were made to any of the experimental counting rates. The 
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FIG. 3: Experimental and theoretical pulse shapes as a function of time, showing EIT, storage and retrieval. 

The color code is: control field - black, pulse in vacuum - blue, delayed, stored and retrieved field - red. Panel (a) with a cw 
control field shows EIT pulse delay. In panel (b) the control field is switched off and then on again after 500 ns, shows light 
storage and retrieval. Panel (c) is similar to (b) but with a 15 fis storage. Panels (d), (e), and (f) are corresponding theoretical 
plots. 

curve fits of g S i are based on a simple theoretical model, and allow us to obtain the efficiency in the idler channel 
and the background contributions to p2 and p% for the stored signal field. These same values are used to produce the 
corresponding theoretical curves in Fig. 4, b and d. The measured values of a < 1, displayed in Fig. 4, b and d, confirm 
the single-photon character of both the source and retrieved signal fields (with the minimum values of a = 0.14 ±0.11 
and a — 0.36±0.11, respectively). Overall, we estimate that the probability p s for successful generation, transmission, 
storage, retrieval, and detection of a signal photon is approximately p s ~ 10 -5 for each trial. The efficiency of photon 
storage and retrieval E can be estimated as the ratio of the values of p2 + pz with and without storage. We find 
E ~ 0.06, in agreement with the theoretical result shown in Fig.3e. 

To investigate the storage capability of our quantum memory at Site B, we measure g S i as a function of the storage 
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FIG. 4: Measured intensity cross-correlation function g s i and anticorrelation function a as a function of the idler 
photoelectric detection probability pi. Panels (a) and (b) are for the source (propagation in vacuum). Panels (c) and (d) 
are for stopped, stored for 500 ns, and retrieved signal field. The solid lines are based on a theoretical model that includes losses 
and background. Error bars represent ± one standard deviation and are based on the statistics of the photoelectric counting 
events. 

10 




FIG. 5: Normalized signal-idler intensity correlation function g S i as a function of the storage time T s at Site B. 

Data (diamonds) are taken for pi = 0.0047, but with a smaller background contribution than that of Fig. 4, c and d. The full 
curve is a fit of the form 1 + Bexp(— t 2 /t 2 ) with B = 7 and the collapse time r = 11 fis as adjustable parameters. Error bars 
represent ± one standard deviation and are based on the statistics of the photoelectric counting events. 



time of the signal field T s (Fig. 5). A Gaussian fit provides a time constant r = 11 fis, which is an estimate of 
our quantum memory time. The collapse is consistent with the Larmor precession of a dark-state polariton in an 
unpolarizcd atomic ensemble in a residual magnetic field |lfil 2lj . Experimentally we attempt to null the uniform, dc 
component of the magnetic field. A definitive way to distinguish whether the collapse is due to uniform or non-uniform 
and ac fields is to measure the damping time of the periodic revivals of the retrieved signal field at longer storage 



times. In a uniform magnetic field, undamped revivals of the dark-state polariton should occur at times equal to nT^, 
where Tl is the Larmor period for level |o) or \b) and n can be either integer or half-integer, depending on the direction 
of the magnetic field relative to the light beam geometry (a synopsis of these ideas is given in the Supplementary 

n 

Information, with the full theory presented in Ref. 29]). We have conducted separate experiments with an externally 
applied magnetic field 30], which suggest that the collapse in the present experiment is likely due to magnetic field 
gradients and/or ac fields at the level of a few tens of mG. However, more extensive investigations to quantitatively 
determine the temporal and spatial structure of the residual magnetic field, and the various contributions to it, are 
ongoing. 

We have demonstrated generation, storage and retrieval of single quanta transmitted between two remote atomic 
ensembles serving as quantum memory elements. The control of the matter-field interface at the level of single quanta, 
and at remote sites, is encouraging for further developments and applications in quantum information science. In 
particular, the storage of a photonic qubit, with two logical states, would represent a crucial advance. In order to 
achieve this, the quantum memory at Site B would likewise need a second logical state, so as to realize a collective 
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If a second logical 



atomic qubit. Two different approaches for such qubits have already been demonstrated 
state were added to both quantum memories at Sites A and B, generation of remote entanglement of two atomic 
qubits would be possible. 
Methods 

To generate single photons at Site A, we use the DLCZ approach in the off-axis, counter-propagating geometry introduced 
by Harris and coworkers 20]. The insets in Fig. 1 indicate schematically the structure of the three atomic levels involved, 
\a), \b) and |c), where {|o); |6)} correspond to the 5Si/2,F — {3, 2} levels of 85 Rb, and |c) represents the {5Pi/2,F = 3} level 
associated with the Di line at 795 nm. The experimental sequence begins with an unpolarized sample of atoms prepared in 
level | a) (left inset of Fig. 1). A 160 ns long write laser pulse tuned to the |o) — ► |c) transition is focused into the MOT with 
a Gaussian waist of about 400 /im. The write pulse generates a cone of forward Raman-scattered signal field via the |c) — > |6) 
transition. We collect a Gaussian mode centered around the momentum k s that forms an angle of about 2° with the write 
beam. The write pulse is so weak that on average less than one photon is scattered into the collected mode for each pulse. The 
signal field is coupled into the 100 meter long fiber connecting Sites A and B. 

For each signal photon emission event, a correlated collective atomic excitation is created in the atomic ensemble. After a 
delay At — 200 ns, a 140 ns long counter- propagating read laser pulse resonant with the |6) — > |c) transition illuminates the 
atomic ensemble and converts the atomic excitation into the idler field. Under the conditions of collective enhancement, the 
idler field is emitted with high probability into the mode determined by the phase-matching condition ki — k w + k r — k s , where 
ki,k w and k r are the wave vectors of the idler, write and read fields, respectively. The waist of the signal-idler mode in the 
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MOT is about 150 /im. The idler field is directed onto a single photon counter Dl. Ideally, photoelectric detection of the idler 
field projects the quantum state of the signal field into a single photon state. The repetition rate of the experiment is 2 • 10 5 
s _1 . Each data point in Fig. 4 involves an average over a time period that varied from several minutes up to 1.5 hours for the 
data point with the lowest value of pi in d. 

To measure the photoelectric detection probabilities pi, P2,£>3, JJ13, P12, JJ23, and P123, the outputs of the detectors are fed to 
three "Stop" inputs of the time-interval analyzer which records the arrival times with a 2 ns time resolution. The electronic 
pulses from the detectors D1,D2,D3 are gated for periods [*o» *o + T g], with Tg = 140 ns, T 9 2 = T 9 3 = 240 ns, respectively, 
centered on the times determined by the write and read (for no storage) or control-retrieval (for storage) laser pulses. Counts 
recorded outside the gating periods are therefore removed from the analysis. The list of recorded events allows us to determine 
the single-channel photoelectric event probabilities = Ni/M, where Ni is the total number of counts in the i-ih channel and 
M is the number of experimental trials, (for Di, i — 1,2,3). If photoelectric detections in different channels i,k,m happen 
within the same gating period, they contribute to the corresponding joint probabilities pij — Nij /M, where Nij is the total 
number of coincidences between Di and Dj, i,j = 1, 2, 3. The joint probability of all three detectors registering a count is given 
by P123 = N 123 /M. 

I. APPENDIX 

Here we provide additional information on the theory of electromagnetically induced transparency (EIT) for atoms 
with Zeeman degeneracy which is essential for adequate modelling of the system. We also include details of the analysis 
of the photoelectron counting statistics of the light fields detected in our experiment, in support of generation, storage, 
and retrieval of single photon states of the electromagnetic field. 

A. Model description 

At Site B we have an atomic ensemble of TV atoms with three hyperfine levels in a A configuration labeled by |o), 
\b), and |c) with total angular momenta F a , Ff,, and F c respectively, as illustrated in Fig. 1. These atoms interact 
with a classical control field on the a — > c transition and a signal field on the b — > c transition. The dynamics of the 
interacting system is governed by the Hamiltonian 

H = H + V AS + Vac + V AB (2) 
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where H = Ek.A^k^kA 6 ^ + E™ c= ^ c ^ ^ I/ 1 ) + fc » E t-F. ^" S W) is the sum of the free 

signal field and atomic energies. The slowly varying operators for atom /x, are defined by o"J* ™* (/i) = 
l s i m i)/i i s 2 m 2\ exp (— i (lu S2 — ui Sl ) (t — Zp/c)). Here Vas = ~ Eu ' -^( r M) ^ s the signal field-atom inter- 
action where D M is the dipole operator for atom fi and E (r„) is the quantized signal field operator at 
the atom's position r M . The atoms are influenced by the classical control field via the interaction Vac — 
""^ (*) Em^-f. 1- m a +r E,! ^ i +r f/ 1 ) + h.c, where fi (i)C^° \ m a +r is tnc ^ ara frequency of the control 
field with polarization r. We also account for the presence of a weak, uniform DC magnetic field through the Zeeman 
interaction 

Vab =VbJ2 {daF^ + fliPjf) + g c F^) ■ B, (3) 

where Ff' = PgF^Ps is the projection of the atomic angular momentum onto the hyperfine level s. As we shall 
demonstrate later, the associated Larmor precession serves as a possible explanation for the reduction of read-out 
efficiency of stored photons over long storage times. 

Before the signal field interacts with the ensemble, each atom is prepared in an unpolarized state. We identify this 
as the mixed atomic "vacuum" state, p va ci as there are no spin coherence excitations present: 

, N F b 

~ E *12»- ( 4 ) 



(2F b + V A 



The signal field is described by the slowly varying photon density operator, $(z,t) — J^p ( r i t) = 
E/3C/3 (i/VVj J2 q ex P (* (q ' r — w c) t) Sq',/3, where q' = q+ (w c /c) z and (3 = ±1 indicates the polarization £±1 = 
=F (1/V2) (x±iy) of the signal field which propagates in the z-direction. The Rabi frequency on the c «-> b transi- 



tion is defined 



by hg = (c D (hj c /2e V) 1 



1/2 



31| . The decay rate of the optical coherence T cb is one half of the 
spontaneous emission rate of the D\ transition, which for 85 Rb is T c = 2T c b w 2tt x 5.75 MHz. It is convenient to 
define several frequencies resulting from the Zeeman shift of the hyperfine levels in a magnetic field, here chosen for 
simplicity to lie in the z-direction: 

S c = 0g c (^-B z ) 
A c6 = {g c -gb) (yBz) 

S a = (P-r)g a ^B z ) 
A ah = (g a - g b )(^^-B z ^j 

where g Q , g b , and g c are the Lande g-factors for the hyperfine levels a, b and c, respectively. 
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B. Field Propagation 



The response of the medium to a time-dependent control held enables storage of the signal pulse within the 
medium. Following the treatment of Fleischhauer and Lukin |2|, but incorporating the Zeeman degeneracy of the 
atomic levels, and assuming a coherent signal held, we hnd that the /3-polarization component of the signal held obeys 
the propagation equation 



d_ 

dz 



®p (A, z) = i*J dA' (s (A - A') + l -^xp (A, A', z)\ $p (z, A') , 



(5) 



where the non-local linear susceptibility is given by 

XP (A,A') = - 



A' + *T. 



<:b 



(6) 



m=-F b 

where d'p (z) is the spatial derivative of the optical thickness dp, which is defined as the negative logarithm of the on 
resonance intensity transmittance of a signal with polarization (3, in the absence of the control held. Explicitly, dp is 
given by the dimensionless quantity 

2 



dp (z) = 6tt/ c6 (j^-J dz'n (z ) — ^ ^ ^ « m+a 



where X 



m 7 p 



c 



F b 1 F c 
m f3 m+fil 



I \ Er. 



a 



F b 1 F c 

m b P m b +P 



D 



I 



D 



(c D aj 



is the 



and f cb = 

fraction of atoms in the excited state that spontaneously decays into the ground level b, and n (z) is the number 
density of atoms. Furthermore, 

1 L mi p (A, A') 



K m ,p = U(A-A')- 



(7) 



/2tz (A + 5c + mA cb + iT cb ) 

is a kernel whose inverse satishes the property / dA" K (A, A")^ 1 (A", A') = / dA" R- 1 (A, A") K (A", A') = 
6 (A — A'). Here A is the frequency space Fourier variable, with A = corresponding to a signal component resonant 
with the c — > b transition. We also have 



L m ,p (A, A') = 



C 



F a 1 -Fc 

m+/3— r r m+P 



2 1 



/ dA"n (A - A") (A „ M(X \ t ^ (A' - A") 
J (A" + {d a + mA ab )) 



(8) 



With this propagation equation, we are able to numerically calculate the propagation of the signal held leading to 
the various features shown in Fig. 3 (we assume equal signal and control polarizations (3 = r in order to have EIT), 
including group delay pulse storage and retrieval. 
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C. Linear susceptibility for constant control field 

In the limit of constant control field, the susceptibility reduces to the form 

, A - |x m ,ff| 2 (A/r e ) 

X/3 (A, 2) re - — (2) ^ 2 • ( 9 ) 

_w - l«/r c | 2 c„^_ F ;, m+/3 -(A/r c )((A/r c ) + i/2) 



Since the Zeeman shifts are small compared to the spontaneous emission rate, they have been ignored in the above 
expression. Notice that the medium will only exhibit EIT with a judicious choice of atomic hyperfme levels and 
field polarizations. If one of the Clebsch-Gordan coefficients multiplying the Rabi frequency f2 vanishes (i.e. there 
is an excited state |c, m + (3) not coupled by the control field to the corresponding ground state \a, m + /3 — r)), then 
a fraction of the atoms will simply absorb the signal light as if there is no control field present. We point out that 
for the level scheme used in the experiment, EIT is achieved only when the signal and control field polarizations are 
equal {ft = r). The solid curves in Fig. 2 are based on Eq.©. 

D. Dark state polariton 



y 



By generalising the analysis of Fleischhauer and Lukin [2| to an unpolarized degenerate three level system, one can 
derive dark-state polariton operators. These operators, which depend on signal and control field polarizations, exist 
only under conditions of EIT. Assuming EIT is possible, for real f2, the dark-state polariton operator of wavevector 
q and polarization (3 is given by 



ssun. 
3 



^«q+(w c /c)z,,3 + i\J 2F\, + l9* J2m C Fa "l \ ~ m+0-r 



<Mq>*) = < 1 . m+ "~ r r m+f = (10) 

\h 2 + j£n\9\ 2 E, 



t+ g-r r 7n + 

^777. fl 777 + 13 

-,j-'„ 1 p; 



where 



Si %+0-r (q) = \J 2 ~^ 1 E eX P H<1 ' r J & a Z+a-r (m) (H) 

ft 

and is the position of atom fj,. The dark-state polariton operator is a linear combination of electric field and a 
specific superposition of hyperfme spin coherences. 

E. Larmor collapse 



In the presence of a magnetic field, a stored atomic collective excitation Larmor precesses out of the dark-state 
polariton mode into orthogonal collective excitations. This causes an apparent decoherence which can be calculated 
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by evaluating the number of dark-state polaritons as a function of time. Assuming that there are exactly N p (0) 
dark-state polaritons initially, and for a magnetic field oriented in the z-direction, we find the number of dark state 
polaritons is given by 



N p (t) = N p (0)Y,T, 



c 



F b 1 F c 



m'+0 



-r r m + ° m ' 



F a 1 F c 
+0- 



i' + 



Er 



°m m+0 
-,F a 1 J.' c 



• cos ((m — m') A a ;,i) 



(12) 



where again A a & = O^B^) (g a — gb). For short storage times, we find 



.Y,,(/) V,,(t)i«. M )( -^r) 2 (A ab t) 2 



where 



v 2 = 



EE 

m m' 



, 1 F c 
m+0 



1 F c 

j rn' + l) 



a + ° m ' + 0- 



n.'+0 



(to — m ) 



(13) 



E, 



For the level scheme used in the experiment (with /3 = r), rj 2 as 2.1. We have used this model to estimate the 
strength of the magnetic field (assuming it is oriented along the z direction), given the time constant r measured 
experimentally in Fig. 5. We find the frequency to be 



Aqfc 

2tt 



14 kHz. 



(14) 



We have used this value in obtaining the theoretical panels in Fig. 3 and find good agreement with the experimental 
observations. Clearly from Eq. Ijl2|) . we can predict the revival of dark-polariton number when A a bt is a multiple of 2ir. 
The revivals occur at t n = 2irn/[(^-B z ) \g a — <#,|] ss tt/[^^-B z ) \g a \]. By contrast, for a magnetic field perpendicular 
to the z axis, we find # = 2vr/ [(^B z ) \g a \] = 2ir/ [(fB z ) \g b \],Le., # = 2*L 

F. Photoelectric counting statistics 

In order to take into account the possibility of sequences of photon pairs, we use a theoretical model based on 
parametric down-conversion, in which the annihilation operators for the idler and signal field are transformed as 

a { ° ut) = cosh(77)aj m) + sinh^fltM, 



fit™ 4 *) = cosh(77)a^ + sinh^a] 



(in) 



(15) 
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FIG. 6: Normalized intensity autocorrelation functions gu (triangles) and g S s (circles for the source, squares for the stored and 
retrieved field) . Uncertainties are based on the statistics of the photon counting events. 

Here r\ is the Raman gain at Site A. We also wish to assess the overall efficiencies and the background levels in our 
experiment. Modelling the background to the signal in terms of a coherent field with average photon number B Sl we 
find that 

_ (l + 2sinh 2 (7y)) + g s ^ 
9sl ~ sinh 2 ^) + B S 

We also determine the anticorrelation parameter a of Grangier et al. 

sinh 2 (r])(4 + 6sinh 2 (?7))+4B s (l + 2sinh 2 (r/)) . . 

a = t . (17) 

(l + 2sinh 2 (77)+ J B s ) 2 

The singles count rates at detectors Dl, D2, and D3 are given by Ri = eisinh 2 (?7), i?2 = |T| 2 e2 sinh 2 ^) and 
i?3 = |i?| 2 e3 sinh 2 (77) (assuming that Ri <C W, where W is the repetition rate of the experiment); T and R are the 
transmission and reflection coefficients of the beamsplitter BS, shown in Fig. 6. 

In the absence of the medium we empirically find negligible background B s . The solid curve in Fig. 4A is based on 
this model, setting B s — 0. We find that the best fit to the data in Fig. 4A is given by ei « 0.039. The solid line in 
Fig. 4B is based on Ea. H17|) with this value of t\. 

For the stored light, we have to account for the fact that a significant fraction of detected signal photons are due to 
background associated with the control-retrieval pulse. By fitting the data of Fig. 4C to Ea. l|ltj|) . we find B s ss 0.08. 
Substituting this value into Ea. 1)17(1 . we obtain the solid curve in Fig. 4D. In order to reduce this background, we have 
performed initial investigations using an optically pumped Rb cell to filter out light at the frequency of the control 
field. In this case we found increased non-classical correlations between the idler and the stored and retrieved signal 
photon, e.g., for T s = 500 ns g S i increased from 8 ± 0.2 to 15.6 ± 1.4. 
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In addition, we measure the intensity autocorrelation functions g ss = Pizl\piPz\ an d gu- These are shown in Fig. 1. 
In order to evaluate the latter, we insert a beamsplitter and additional detector D a into the path of the idler photon, 
so that gu = p\ a /\piPa[- Using these together with the measured values of g S i shown in Fig. 4 of Rcf.[? ], one can 
evaluate Clauser's parameter R = g1i/[g S s9n\- For classical fields R < 1, whereas we observe strong violation of this 
inequality. 

The total measured transmission and detection efficiencies for the idler and signal fields respectively are Wi = 
0.25 ± 0.03 and w s = 0.15 ± 0.02, consisting of the quantum efficiencies of the detectors 0.55 ± 0.05 and the passive 
transmission losses accounting for the rest. The ratio of h = ei/wi = 0.16 indicates the strength of the spatial 
signal- idler correlations in our source of conditional single photons at Site A, with h — > 1 for the ideal case. 
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